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SIFASP 



A General Computer Program for Simultaneous Factor 
Analysis in Several Populations 

1. Introduction 



1.1 The General Model 

We shall describe a computer program for simultaneously factor analyzing 
dispersion matrices obtained from independent groups. A common situation, 

■when this program ■will be useful, is when a battery of tests has been 
administered to samples of examinees from several populations and one wants 
to study similarities and differences in factor structures between the dif- 
ferent populations - The most important feature of the program is that param- 
eters in the factor analysis models (factor loadings, factor variances, fac- 
tor covariances, and unique variances) for the different populations may be 
assumed to be known a priori or specified to be invariant over populations. 
Given such a specification, the model is estimated by the maximum likelihood 

p 

method yielding a large sample X test of goodness of fit. By computing 
several solutions under different specifications one can test various hypothe- 
ses- For example one can test the hypothesis of an invariant factor pattern. 
The method is capable of dealing with any degree of invariance, from the one 
extreme, where nothing is invariant, to the other extreme, where everything 
is invariant. A detailed account of the method, on which the program is 
based, is given by Jtireskog (1970). 

Consider a set of m populations. These may be different nations, or 
culturally different groups, groups of individuals selected on the basis of 

'^Research reported in this paper has been supported by grant NSF-GB-12959 
from National Science Foundation. 
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some known or unknown selection variable, groups receiving different treat- 
ments, e^c. In fact, they may be any set of exclusive groups cf individuals 
that are clearly defined. It is assumed that a battery of p tests has 
been administered to a sample of individuals from each population. The 
battery of tests need not be the same for eh group, but to be interesting, 
it is necessary that some of the tests in each battery are the same or at 
least content -wise equivalent. 

Let x be a vector of order p , representing the measurements ob- 
g 

tained in group g ♦ We regard x as a random vector with mean vector 

g 

u and variance -covariance Z . It is assumed that a factor analysis 
g g 

model holds in each population so that x can be accounted for by k 

g 

common factors f and p unique factors z , as 
g g 



(i) 



x = [X + A f 
g g g g 



g 



with e(f ) = 0 and S(z ) =0 and A a factor pattern of order p^ x k 
g g g g £ 

The usual factor analytic assumptions then imply that 



( 2 ) 



S = A $ A' + * 2 
g g g g g 



where $ is the variance-covariance matrix of f and \|r is the diagonal 
£ g g 

variance -covariance matrix of z 

g 

In addition to assuming that a factor analytic model holds in each 

population the model may specify that certain parameters in A , $ , \|r , 

g g g 

g = 1, 2, ...,m have assigned values and that some set of unknown elements in 

A , § and \lr are the same for all g . Thus, parameters in A , $ 

g g g g g 

and i|r , g ~ 1,2, ...,m are of three kinds: (i) fixed parameters which 

g 
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have "been assigned given values, (ii) constrained parameters 'Which are unknown 
"but equal to one or more other parameters and (iii) free parameters which are 
unknown and not constrained to "be equal to any other parameter. Equality 
constraints "between parameters for the same populations may also "be used 
though this would he unusual in practice. The advantage of this approach 
is the great generality and flexibility obtained by the various specifica- 
tions that may be imposed. The most common situation is when the same 
battery has been administered to each group and when the whole factor 

pattern A is assumed to be invariant over groups. This case will 
§ 

hereafter be referred to as the standard case* 



1,2 Identification of Parameters 

Before an attempt is made to estimate a model of this kind, the iden- 
tification problem must be examined. The identification problem depends 
on the specification of fixed, free and constrained parameters. Under a 
given specification, reach A , $ and If generates one and ^nly one 



Z but it is well kno'^ th' 1 * 1 
g 



1 A and © generate the sen 



Z . It should be noted that if A is replaced by A T and $ by 
g g g g g 

T $ T T , where T is an arbitrary nonsingular matrix of order k k , 
g g g g B I 

2 

then Z is unchanged. Since T has k independent elements, t iis 



suggests that 



2 



independent conditions should be imposed on A and/or 



m 2 

$ to make these uniquely aefined and hence that Z k independent 
g g=i g 

conditions altogether should be imposed. However, when equality constraints 

over groups are taken into account, all the elements of ell the tran forma- 
tion matrices are not independent of each other and there "'ore a lessor 
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nuiriber of conditions need to be imposed. it is hard to give further spe- 
cific rules in the general case. To make sure that all indeterminacies 
have been eliminated, one should verify that the only transformations 

that preserve the specification about fixed, free and con- 
strained parameters are identity matrices. 

In the standard case when the whole factor pattern is invariant over 
groups, however, a more precise consideration of the identification problem 
can be given. Suppose that the A is replaced by A* » AT* 1 and each § 

g 

is replaced by # = T^> T T , g = 1,2, ...,m , where T is an arbitrary 
§ g 

nonsingular matrix of order k x k . Then each Z remains the same. 



g 

2 

Since the matrix T has k independent elements, this means that at 
2 

least k independent conditions must be imposed on the parameters in 

A , 10 ma ^ e ‘kttese uniquely defined. 

The most convenient way of doing this is to let all the $ be free 

g 

and t- fix one nonzero olement and at least k - 1 zeros in each column 
of A . In an exploratory study one can fix exactly k - 1 zeros in 
almost arbitrary positions. For example one may choose zero loadings 
where one thinks there should be "small" loadings in the factor pattern. 
The resulting solution may be rotated further, if desired, to facilitate 
better interpretation. In a confirmatory study, on the other hand, the 
positions of the fixed zeros, which often exceed k - 1 in each column, 
are given a priori by an hypothesis and the resulting solution cannot be 
rotated without destroying the fixed zeros. 
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1.3 E s t imat ion and Testing of the Model 

th 

Let N be the number of individuals in the sample from the g popu- 

lation and let x be the usual sample mean vector and S g the usual sample 

variance -covariance matrix with n = N - 1 degrees of freedom- The only 

& & 

requirement for the sampling procedure is that it produces independent 

measurements for the different groups- 

If we assume that x has a multinormal distribution it follows that 

§ 

S has a Wishart distribution based on Z and n degrees of freedom, 
g g s 

The logarithm of the likelihood for the sample is 

(3) log L g = ~§ n g [log|2 g ! + tr(S g S' g )] - 

Since the samples are independent, the log-likelihood for all the samples 

m 

(1^) log L - £ log L 

S=1 

Maximum likelihood estimates of the. unknown elements in ; ^g 9 

g = 1^2, ... ; m , may be obtained by maximizing log L • However, it is 
slightly more convenient to minimize 

(3) F = § z^ n g [log|2 g | + ^(Sg^g 1 ) - losing! - P] 

instead. At the minimum, F equals minus the logarithm of the likelihood 
ratio for testing the hypothesis implied by the model against the general 

alternative that each Z is unconstrained- Therefore, twice the minimum 

® 2 
value of F is approximately distributed, in large samples, as "X- with 

degrees of freedom equal to 




S 
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d = | ?(P + 1) - t 

where t is the total number of independent parameters estimated in the 
model. 

The minimization of F with respect to the independent parameters is 
done by means of a modification of the iterative method of Fletcher and Powell 
(1965) described by Gruvaeus and j£3reskog (l97°)* The minimization method 
makes use of the first-order derivatives and approximations to the second- 
order derivatives of F and converges rapidly from an arbitrary starting 
point to a local minimum of F - If there are several minima of F there 
is no guarantee that the method will converge to the absolute minimum. 

The adaptation of the problem of minimizing F to the Fletcher- 
Powell method is described by Jtireskog (1970; section 2.4). 



1,4 Scaling of Observed Variables 

VJhen the units of measurements in the different tests are arbitrary, 
it is usually convenient, though not necessary, to rescale the observed 
variables, before the factor analysis. Let 



m 



(T) 


S = 


(l/n) S n £ 
g=l S 




m 




with 


n = L 


n and let 




g-1 


g 


(8) 


D = 


(diag S) 1 ^ 2 



Then the variance -covariance matrices for the rescaled variables are 
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( 9 ) S* = BS D , g = 1,2, ...,m • 

& 

The weighted average of the S* is a correlation matrix* The advantage of 

§ 

this rescaling is that, when combined with the rescaling of the factors as de- 
scribed in the next section, the factor loadings are of the same order of magni- 
tude as usual when correlation matrices are analyzed and when factors are stan- 
dardized to unit variances- This makes it easier to choose start values for 
the minimization (see Jiireskog, 1970, section 3*5) and interpret the results - 
It should be pointed out that it is not permissible to standardize the 
variables in each group and to analyze the correlation matrices instead of 
the variance -covariance matrices- This violates the likelihood f unct i on (k ) 
which is based on the distribution of the observed variances and covariances. 

1.5 Scaling of Factors 

The fixed nonzero loading in each column of A can have any value- 
This is only used to fix a scale for each factor which is common to all 
groups- In the standard case, when the maximum likelihood solution has 
been obtained, the factors may be rescaled so that their average variance 
is unity. This rescaling is obtained as follows. Let 

m 

( 10 ) $ = (l/n) 2 n 8 , 

g=l 2 S 

m 

with n = 2 n .as before, and 

g=l « 

(11) D = (diag $) 1 / 2 
Then the rescaled solution is 
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. x /N ^ -1 

(12) A* = AD 

(15) 3* = D$ D , g = l,2,...,m . 

§ § 

s\ y\ 

The matrix A* has zeros wherever A has zeros "but the fixed nonzeros 

A A 

in A have changed their values. The weighted average cf ^ is a 

8 

correlation matrix. 
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2. The Program 

In this section we describe briefly what the program does- Details 
about the input and output are given in sections 3 and k respective! y. 



2*1 What the Program Does 

The input data may be correlation matrices with standard deviations 
or dispersion matrices. From these input matrices, variables may be 
selected to be included in the analysis, so that the matrices to be 
analyzed may be of smaller order than the input matrices- Variables may 
also be interchanged with one another. The matrices to be analyzed may be 
dispersion matrices or dispersion matrices scaled by the program (see 1.4). 

The user can request an accurate or an approximate solution. If an 

accurate solution is requested, the iterations of the minimization method 

are continued until the minimum of the function is found, the convergence 

criterion being that the magnitude of all derivatives be less than . 00005 N, 

m 



where 



N = 



(l/m) 



Z n 
g =i s 



The solution is then usually correct to three 



significant digits. If an approximate solution is requested, the iterations 

terminate when the decrease in function values is less than 3^* The approxi- 

2 

mate solution may be useless but the residuals and the value of X will 
usually give an indication of how reasonable the hypothesized model is- The 
option of an approximate solution has been included in the program for the 
purpose of saving computer time in exploratory studies where the primary 
purpose is to find a reasonable model. Once such a model has been found^ 
an accurate solution may be computed. 
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A variety of options for the printed output is available. Residuals 

t for each population may be printed* These are defined as the differences 

between observed (S ) and estimated (S ) variances and covariances, which 

§ § 

* 2 

are useful for judging the goodness of fit of the model to the data. X is 

printed as an overall goodness of fit test statistic and, in one version of 

the program, standard errors for the estimated parameters may be requested 

(see 2*3)« 



2.2 How Fixed , Free and Constrained Parameters Are Specified 

The elements of the parameter matrices are ordered as follows. The 
matrices are assumed to be in the order A^,Ag, ... ; A , ^2* * * * 9 ^m 9 

^1^ ^ 2 9 * * * 9 an( ^ eacl1 matrix, the elements are ordered row-wise. 

Only the lower half including the diagonal of the symmetric matrices 

are stored. The diagonal matrices • . . .\lr are treated 

as row-vectors. 

for each of the parameter matrices, a pattern matrix is defined, with 
elements 0, 1, 2 and 3 depending on whether the corresponding element in 
the parameter matrix is fixed, free, constrained follower and constrained 
leader, respectively. A constrained parameter is called a constrained 
leader the first time it .appears in the sequence. The parameters, appearing 
later in the sequence and assumed to be equal to the constrained leader are 
called constrained followers. 

The above technique defines uniquely the positions of the fixed, free 
and constrained leader parameters. It does not define, however, which 
followers go with which leaders, if there is more than one leader. To do 
so one must specify all the followers associated with a given leader. This 
is done by assigning to each leader and follower a four-digit number MCCC, 

O 
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where M defines the matrix in which the constrained parameter appears. 

M = 1 for A , 2 for <2> and 3 for y where A is A^, A^ . . . ^ A^ reading 
row-wise one matrix after the other^ $ is 0 ...,$ and ^ is 

V n ytypj • - *■ • The position of the parameter in the matrix is described by 

CCC . For example ^ 

1001 1005 2003 

defines the first element in A , ; to be equal to the fifth element in 



A A as well as the thira elemeno zj 

5 



;> 9 



where A. is the 



leader and A_ and 

5 



are the followers. 

3 



Pattern matrices have to be provided for each matrix containing both 
fixed and free parameters and for each matrix containing constrained param- 
eters- Patterns for matrices whose elements are all fixed or all free are 
set up by the program. 

We give a simple example to illustrate the above specifications. 
Suppose we have two populations and 



A = 
1 



A 


0 








f. 


0 


0 


0 " 


1 










1 








A 


0 


fc = 


1 


\|f = 


0 


'dr 

*2 


0 


0 


0 


A 6 


1 


l 

r— 1 

CM 


*1 


0 


0 


i dr 
3 


0 


0 


A 8_ 








0 


0 


0 


1 

** 



A = 



K 


0 








f 


0 


0 


0 


9 










5 






A 


0 




"l 




0 


f 


0 


0 


11 




II 

CM 

O 




*2 = 


6 


0 


A li, 


.*5 \ 


0 


0 


'dr 

7 


o ! 


0 


A i6_ 








0 


0 


0 


\ 
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and the specifications of leaders and followers are 



1001 


1003 


1009 


1011 


10 Ik 


1016 






3001 


3002 


3005 


3006 


3007 


3008 







In this model ten independent parameters will be estimated. This is the number 

of. and 1 ! s in the pattern matrices. 

In addition to the above specifications for fixed, free and constrained 

parameters, start values have to be given for all parameters, except when 

one or more of the parameter matrices are of standard form, i.e., A = I , 

$ = I ? ^ = 0 , g = l, 2 ,...,m. The start values define the fixed param- 
& & 

eters and initial values for the minimization procedure for the other param- 
eters. Constrained parameters which are assumed to be equal must be given 
the same values. Otherwise, initial values may be chosen arbitrarily but 
the closer they are to the final solution the less computer time it will 
take to reach this solution (see Jtireskog, 1970, section 5-5). 



"ith \ " S - S " V ' *14 = A iS ' *l - *2 - *5 - *6 “ d *T 

The pattern matrices for , ® 1 , 1> 2 , ^ and \|r are 






3 o 
2 0 
0 1 
0 i 



u 



P = [3 

+1 



P. = 



^2 



r 2 o' 

2 0 
0 3 
0 2 



P 4> = 

2 



* [°1 o ] \ - [2 2 3 



O 
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2.3 Limitations 

The program is written in PORTRAIT TV-G and has been tested out on the 
IBM 360/65 at Educational Testing Service. Double precision n. used in 
floating-point arithmetic throughout the entire pregram. With minor changes 
the program should run on any computer with a FORTRAN IV compu r. In 
computers with a single word length of 36 bits or more, single precision 
is probably sufficient. 

Three versions of the program are available: SIFASP, SFASPL and SFASPF. 

Their limitations as to the maximum number of populations,, variables, factors 
and independent and nonfixed parameters they can handle as well as their 
storage requirements on the IBM 360/65 are given in the following table* 

The given storage requirements assume the programs are overlayed. 







SIFASP 


S FAS PL 


SFASPF 


Max* no. of populations (m) 


10 


10 


10 


Max. no. of 


variables (p) before selection 


120 


200 


120 


Max* no. of 


variables (p) after selection 


24 


40 


24 


Max. no. of 


factors (k) 


12 


20 


12 


Max. (p-p(p 


+ 1)) 


312 


820 


312 


Max. (mpk) 




288 


800 


288 


Max. (mp) 




48 


80 


48 


Max. (— k(k 


+ 1)) 


78 


210 


78 


Max. no. of 


independent parameters 


120 


200 


120 


Max. no. of 


nonfixed parameters 


150 


300 


120 


Storage requirements (K = 1024 bytes) 


i 44 k 


280K 


i 46 k 



SlFASP and SFASPL.&re identical except for dimensions. Neither of 
these programs use expressions for second-order derivatives; instead the 
matrix of the Fletcher and Powell procedure is an identity matrix 

(see Fletcher & Powell, 1963; Jttreskog, 1970; or Gruvaeus 8 c jBreskog, 
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1970)- SFASPF, on the other hand, makes use of such expressions and the 
speed of convergence is therefore somewhat faster* Standard errors for 
the estimated parameters can only be obtained with SFASFF. 

2*4 Availability 

A copy of the program may be obtained by writing to one of the authors* 
The user must provide a tape on which the program will be loaded* The program 
will be written on the tape with 80 characters per record* The tape will be 
unlabeled* The user must specify whether he wants the tape blocked or un- 
blocked, on 7 -track or 9-track, in EBCDIC or BCD mode, as well as the density 
and parity required* Test data will be at the end of the program* The test 
data are described in the Appendix. Anyone using the program for the first 
time should make sure that the test data run correctly. 

2*9 Disclaimer 

Although the program has been working satisfactorily for all data 
analyzed so far, no claim is made that it is free of error and no warranty 
is given as to the accuracy and functioning of the program. 
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3* Input Data 

For each data to he analyzed, the input consists of the following •» 

1- Title card 
2 Parameter cards (2} 

3* Selection of variables from the input matrix 
km Input matrices 

5* Pattern matrices for the parameter matrices 
6. Equalities 

7* Initial values for the parameter matrices 
8. New data set or a STOP card 

Sections 3*1 through 3*8 describe in general terms the function and setup 
of each of the above quantities* Illustrative examples are given in the 
Appendix* 

"Whenever a matrix or vector for m populations is read in it is pre- 
ceded by a format card, containing at most 80 columns, beginning with a left 
parenthesis and ending with a right parenthesis. The format must specify 
floating point numbers for the input and parameter matrices, and fixed 
point nu mb ers for the pattern matrices, consistent with the way in which 
the elements of the matrix axe punched on the following cards. Users who 
are unfamiliar with FORTRAN are referred to a FORTRAN Manual, where format 
rules are given. Matrices are punched as one long vector, reading row - 
wise, each population beginning on a new card . For the symmetric matrices 
only the lower half of the matrix including the diagonal should be punched . 
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5-1 Title Card 

■Whatever appears on this card will appear on the first page of the 
printed output. All 80 columns of the card are available to the user. 



cols. 6-10 

cols. 11-15 
cols. 16-25 



5.2 Parameter Cards (2) 

Card 1 : All quantities on this card; except for the logical indicators, 

must be punched as integers right adjusted within the field, 
cols. 1-5 Number of populations m 

Order of the input matrix ( p ), before selection of 
variables 

Number of columns in A (k) 

Total estimated execution time in seconds for all 
stacked data (SEC). This should be a number slightly 
less than the time requested on the control cards so 
the program will have time to print and/or punch results 
up to that point. (Note: SEC should be read in for 

each data set and should be the same for all data sets 
in the stack.) 

Logical indicators (see below) 

Integer output indicators (see below) 

Logical Indicators (cols. 51-37): The logical indicators control the 

input and output as described below. 

Column 51 determines whether dispersion matrices, or correlation matrices 
and vectors of standard deviations, are read in as input to determine the 
matrices to be analyzed. 



cols. 31-37 
cols. ^5 -^*6 
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col* 31: r= T , if a dispersion matrix -with diagonal is read in for 

each population 

col* 31: = F , if correlation matrices without diagonal^ followed by 

vectors of standard deviations are read in for each 
population 

Column 32 determines whether the matrices , g = l>2,...,m to be 

— ■ & 

analyzed are different from the matrices analyzed in the previous data set. 
col* 32: = T , if new matrices are to be analyzed (note: this is 

always true for the first data set) 
col- 32; = F , same matrices as for previous data set are analyzed 

Column 33 determines whether the matrices to be analyzed are scaled 
or not * 



col* 33: » T , 



matrices to be analyzed are scaled by the program to 

S* = DS D , g = 1,2, • where D = (diag S) , 

g § 



S 




S 

g 



n 



m 

2 

g=l 



n 

g 



col. 33: « F , analysis performed on the unsealed S , g - 1,2, •*.,m 
Column 34 determines whether selection of variables from the input 
matrices is desired. 

col* 3^: « T , if selection of variables is wanted 

col* 3^: = F , if no selection of variables is wanted 

Column 33 determines whether we are considering the standard case or 



not • 

col. 35: = T , the standard case is considered (i.e., A , i = 

are constrained to be equal to A- ; in this case the 
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pattern matrix and starting matrix for t. will be 
read in for the first population only) 
col. 55: = F , we are not considering the standard case (pattern 

matrices and starting values for all m populations 

■will he read in) 

Column 56 determines whether the starting values for , g - 1,2, ...,m 

are dispersion matrices or correlation matrices with standard deviations from 

which the dispersion matrices will he computed. 

col. 5 6: = T , starting $'s are dispersion matrices 

col. 56: = F , starting 4’s are correlation matrices without diagonal 
a-.-iri with standard deviations 

Column 7T determines whether an accurate or an approximate solution 
is required* 

co ^ jY* s T , if an approximate solution is required 
col. 57: = F , if an accurate solution is required 

Integer Output Indicators (cols- 1+5-46) 

Column 1+5 determines the type of printed output wanted. This can he 
standard output ( S ), the matrices to he analyzed and parameter specif nca 
tions ( K ), residuals and S for each population ( C ), and technical output 

from minimization ( T )• 

col. 1+5: = 0 ; for S 



col. 1+5 
col. 45 
col. 45 



= 1 , for S + R 

= 2 , for S + C 

= 3 } for S + R + C 



col. 45: = 4 j 
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col. 45! = 5 , for S + R + T 

col. 4-5! = 6 , for S + C + T 

col. 45: = 7 , for S + R + C + T 

Column 46 determines certain extra printed or punched output. This 
can he standard errors ( P ) which is only applicable to SFASPF, punched 
solution ( P ), and a scaled solution ( G ) 



col. 


46: 


= 


0 


> 


if no extra output 


is 


•wanted 


col. 


46: 


= 


1 


> 


for 


F 


(never set 


to 


1 for SIFASP or SFASPL) 


col. 


46: 


= 


2 


> 


for 


P 








col. 


46: 


= 


3 


> 


for 


F 


+ P 






col. 


46: 


= 


4 


> 


for 


G 








col. 


46: 


= 


5 




for 


F 


+ G 






col. 


46: 


= 


6 


? 


for 


P 


+ G 






col. 


46: 


= 


7 


; 


for 


F 


+ P + G 






Card 


2: 


This 


card -will 


specify the 


number of observations or sample 



size for each population. Thus there will be m integer numbers punched, 
right -adjusted in five column fields. 

Caution : Mien specifying m , p , h on card 1 of the >■ rameter 

cards be sure you have read the limitations imposed on them (see 2.3)* 

3. 3 Selection of Variables 

These cards will be read in only if the parameter card has a T in 

column 32 and a T in column 34. Omit otherwise. 

The first card will have an integer value P new punched in columns 

1-5, right adjusted within the field. This integer will specify the order 

of the S , g =. 1,2, ...,m after selection (P new S P) * 

§ 
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The next card will contain integers, right-adjusted in five column 
fields, (i.e., sixteen such values will fit on one card) specifying which 
columns (rows) are to be included * For example: if P = ° > P n evj “ ^ 
and the 1st, 2nd and 5th columns (rows) are to be excluded . This card 
would have a 3 punched in column 5; a ^ punched in column 10 and a 6 
punche . in column 15- 

Note that if p = p there will be no reduction in the size of 
new ■ 

the S but columns (rows) can be interchanged. 

Input Matrices 

Omit if column 32 of the parameter card is F . Otherwise read in a 
format card followed, on subsequent cards by the input matrices, starting a 
new card for each population. 

If column 31 of the parameter card is F , the input matrix for the 
first population, preceded by a format card, is read in without the diagonal. 
This is immediately followed by a format card and the vector of standard 
deviations for the first population. Subsequent cards are input matrices 
without diagonal for the remaining populations each followed on a new card 
by its vector of standard deviations, and starting a new card for each 
population. The formats for the first population will apply to subsequent 
populations . 

3*5 Pattern Matrices 

The pattern matrices are preceded by a data card with entries in columns 
1-3, the column defining the matrix in question, 1 for A } 2 for 0 and 
3 for \|r . 
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cols. 1-3: CCC where C = 0 , if the matrix is fixed 



C = 1 , if the matrix is free 



0 = 2 . if the matrix has mixed values 



A pattern matrix should be provided only when C = 2 (see 2.2). 

For example, if columns 1-3 are punched 20 1, the matrix A (i.e.. 



The pattern matrix consists of a format card specifying an I -format 
and subsequent cards with the integer entries of the parameter matrix, 
beginning a new card for each population. 

3.6 Equalities 

Omit if the pattern matrices do not contain any elements 2 or 3« Other- 
wise starting in column 1 punch the four-digit numbers MCCC as described 
in section 2.2. For each new constrained leader start a new card. The last 
entry on each "equality” card is a zero indicating more "equality" cards 
follow, or a four indicating it is the last one (see Appendix for examples). 

3.7 Initial Values for the Parameter Matrices 

The initial values are preceded by a data card with entries in columns 
1-3, the column defining the matrix in question. 

cols. 1-3: CCC where C = 0 , if the matrix is of standard form 




1, 2, ...,m are read in* 



® 4> g , g = 1,2, ...,m ) 



, m ) is all free. In this case 



( see 2.2 ) 



C = 1 , if the matrix is nonstandard 



This card is then followed by the necessary start values (see 2.2) for matrices 



with C = 1 • 



That is, each nonstandard matrix of m populations is read in 
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with its own format card^ starting a new card Tor each population. If column 
36 of* the parameter card is F , , preceded by a format card, is read in 

without the diagonal. This is immediately followed by a format card and the 
vector of standard deviations for the first population. Subsequent cards are 
> i = 2, 3j • • without diagonal each followed on a new card by its vec- 
tor of standard deviations, and starting a new card for each population. The 
formats for the first population will apply to subsequent populations. 

3.8 Stacked Data 

In sections 3.1 to 3 • 7 we have described how each set of data should 
be set up. Any number of such sets of data may be stacked together and 
analyzed in one run. After the last set of data in the stack, there must 
be a card with the word STOP punched in columns 1-^. 
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4. Printed and Punched Output 

The output consists of a series of printed and punched tables as described 
in section 4.1-4. 7. Examples of printed output are given in the Appendix. 

4.1 Standard Output (s) 

The standard output is always obtained, regardless of the value punched 
in columns 4-5 and 46 of the parameter card (see 3*2)* The standard output 
consists of the title with parameter listings the final solution and the 
result of the test of goodness of fit. 

The parameter listing gives the information supplied on the parameter 

card. 

The final solution consists of the three matrices A 9 $ and \|r , 

printed for each population. 

2 

The test of goodness of fit gives the value of X and the corresponding 

degrees of freedom. The probability level is also given. This is defined 

2 

as the probability of getting a X value larger than that actually obtained, 
given that the hypothesised structure is true. 

Just above the table giving the final solution, the following message is 
printed 



"im = X" 

Usually X - 0 , but if, for some reason, it has not been possible to determine 
the final solution., X will be 1, 2, 3, 4 or 5 . If XND is 1, 2 or 3* "seri- 
ous problems" have been encountered and the minimisation of the function 
cannot continue. One reason for this may be erroneous input data. Another 
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reason may be that a point has been founds -where one of the matrices £ is 

§ 

not positive definite* A third reason may be that insufficient arithmetic pre- 
cision is used. If IND is 4, the number of iterations has exceeded 250. If 

IHD is 5/ the time limit SEC has been exceeded (see 5-2). If HID / 0 , 

the solution obtained so far is automatically punched on cards in such a 

■way as to be Immediately available as initial estimates for a new run with 

the same data. Thus there is little loss of information when execution is 
terminated with HID / 0 * 



4.2 Matrices S and Parameter Specifications (r) 



If column 45 of the parameter card is 1, 3 ; 5 or 7 (- " ~-2), the. 

matrices to be analyzed;, S , g = 1^2, ••• > m . as obtained! after exclusion 

§ 

of variables and/or scaling (see 1.4),, if any., are printed.- These matrices 
are printed row-wise with four decimals* A1 o a table of parameter speci- 
fications^ containing the information provided by the pattern matrices (see 
2.2), is printed. For each population, three integer matrices are printed 
corresponding to A 9 $ and \|r . In each matrix an element is an integer 

equal to the index of the corresponding parameter in the sequence of inde- 
pendent parameters. Elements corresponding to fixed parameters are 0 
and elements corresponding to the same constrained parameter have the same 
value. Examples are given in the Appendix. 



4p Technical Output (T) 

If column 45 of the parameter card is k, 5, 6 or 7 (see 3.2)., the 
technical output is printed. This consistr of a series of tables which 
describe the behavior of the iterative procedure and give various measures 



-25- 



of the accuracy of the final solution- Ordinary users will have little 
interest in these tables. 

The first table of the technical output gives the initial estimates for 
Ag > > ^g > g = 1 , 2 , „ ..,m . 

The next two tables show the behavior of the iterative procedure under 
the steepest descent iterations and under the following iterations by the 
Fletcher and Powell method. For interpretation of these tables the reader 
is referred to Gruvaeus and JBreskog ( 1970 ). If something gees wrong, so 
that TNT) is 1, 2 or 3 (see 4.l), these tables may contain valuable 
information • 



4 A Matrices £ and Residuals 

S 



(C) 



If column 45 of the parameter card is 2, 3^ 6 or 7 (see 3*2), the 

A a ^ a ^2 A 

matrices £ = A $ A* + and the residual matrices S - £ , g = 1, 2, . ..,m , 

o g g g g g § 

/\ 

are printed. The matrices £ are computed from the final solution* If 

g 

the fit is good, £ should agree well with S and the residual matrices 

g g 

should be small. Elements of the residual matrices may suggest how the 
hypothesized structure should be modified to obtain a better fit. The 
matrices are printed row -wise, each element with four decimals. 



4.5 Scaled Solution (G) 

It column 46 of the parameter card is 4, 6 or 7 (see 3.2), a scaled 

solution is printed. (See 1*5 on scaling of factors.) 
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i+»6 Standard Errors (F) 

If column U-6 of the parameter card is 1, 3 ; 5 °r 7 (see 3-2). large 

sample approximations to the standard errors of the estimated parameters 
are printed. These are printed row-wise in matrix form and each number 
is printed -with three decimals. The reader is referred to the paper by 
jSreskog (1970 ) for infcmnation about ho*w the standard errors are obtained. 

The standard errors are for the parameters of the unsealed solution. 

Punched Output (P) 

If column k6 of the parameter card is -2. p, 6 or 7 ( see 3*2), the 
final solution is punched on cards • The matrices are punched on cards in 
rector form, reading row -wise, beginning a new card for each population 
and each matrix. Each of the three matrices A , $ , are preceded 

by a format card where by A we mean A , g = 1, 2, ...,m , by 0 we mean 

o 

$ , g = 1,2, . ..,m and by \[r we mean , g = 1, 2, • . «.,m • In the 

standard case only one A is punched regardless of the number of populations* 
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APPENDIX 



; shall -ill ir-Rtra- - i^P’U 13 data are set up and what the printout 

looks —ike by means of Uhc? sf^all sets of* data. These data also serve as test 
data to he run when the ^tfc-gTsun has been compiled on a different computer. 

Both sets of data analyzed in one run with SPASFF- Pages A*+-A5 

show card by card hc'“ input data are punched. One line corresponds to 

one cax*.i. Pages A6-Af.‘l sho^ r the corresponding printout obtained. 

The first set of "Holzinger -Swineford Data/' consists of four 

9x9 correlation matrices without diagonal each with a set of standard de- 
viations. All variable^ #09e Included in the analysis^ and the input matrices 
are to be scaled by the before being analyzed. The following model 

is assumed: 



(9x5) = Ag *= *= — 



*5 + 6 . 
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(V 

5 Ip 

! 4) cb ^ 

L -7 18J 



0 >,. 




11 

£3 *2b 



and the \|r* s are constrained to be equals 



« V 2 - V-, = % = (♦ 1 ,VVW^V t 8 ,t 9 ) * 

Initial values for A , $ and \|r were obtained from preliminary analyses 

S g S 

of each population separately. All printed output is requested. 

The second set o_" data, "Artificial Data for Illustrative Purposes, " 

th 

consists of two 10 x 10 dispersion matrices with the 10 variable excluded 
and with the vaxxable moved to the first position. Ine following model 

is assumed: 
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<t> 




1 




7 


*1 - 


*2 h 


<3> = 

" 2 


*8 <t> 9 




<t>. <t> 

456 




ho hi 0 i2_ 



U = [ f 2 , f y f 6 , f Q , t g ] , 

h = ^10^11^12^13" ^14'' ^15^16^17^18 •* 

In this analysis we impose the constraints 7\^ = 7\^ ^ T^q = 7 \^ , = 7^] 

7 \^q = "A^ • Initial values have been chosen as 0.9 ^ or a-11 nonfixed Vs 

except 7\g and which have an initial value of 0.4. All vjr T s have 

initial values of 0.8. The $ are read in as correlation matrices, each 

6 

followed by its vector of standard deviations. Only the standard output, 
the matrices to be analyzed and the parameter specifications and the standard 
errors are requested as printed output . 

At various places in the output, time estimates are printed. The time 
shown is the time taken to compute the solution that follows the time esti- 
mate. This time includes only the iterations and not the time for printing, 
except possibly the technical printout. 
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0.274 


6. 280 


0.7 84 



POPUL ATI O N 3 



S 



12 3 A 5 6 7 6 9 

1 0.923 

z 0 . 199 0 . 745 

3 0.219 0.180 0.982 

A 0.297 0.042 0.220 0.931 

3 0. 357 O. 211 O. 1R9 0.696 1. 12 1 _ 

6 0.349 0.007 0.110 0.575 0.70? 1.0l8 

7 0.214 0*009 -0.070 0.08B 0.118 0.123 1.027 

8 -0.021 -0.010 - 0*143 0.054 0.094 0.034 Q..2J.4 1.230 

Q 0*086 -0.121 -0.059 0.154 0.021 0.121 0.152 0.321 C^998 



P OPUL A T I_ON_. 



S 



1 2 3 4 5 6 7 8 9 

1 1.126 

2 0.410 1. 46 1 

3 0. 523 0.409 1.053 

4 0.308 0.013 0.064 1.072 

5 0.293 0. 01 3 0. Oil 0.824 1 .12 7 _ 

6' 0.420 ' 0.311 0. 102 0.615 0.662 0*980 

7 0.523 0.454 0.265 0.182 0.087 0.418 1.377 

9 0.128 0. 061 Q. 031 -0.033 0*064 0*189 0^3 4,1 1.0 04 

” 9 "" 0.255 - 0.050 0.011 -0.054 O.lll 0.248 0.233 0.398 1.091 




*1 r™ 



PARAMETER SPECIFICATIONS 
POPULATION 1 



LAMBDA 

0 0 0 - 

10 0 
2 0 0 

0 O . _ - 0 - 

0 3 0 

0 4 0 

o o a 

0 0 5 

0 0 6 



PHI 

7 

8 9 

10 11 12 

PS1 . 

13 14 15 16 17 IP 19 20 21 

POPULATION 2 

L AM BOA 

_ .0 -Q -0 - 

1 0 0 

2 0 0 

POO 

0 3 0 

0 4 0 

0 0 Q 

0 0 5 

0 0 6 



PHI 

22 

23 24 

25 26 27 

PST. - 

13 14 15 16 17 18 19 20 2l 

P OPULATION 3 



lASiirtft 

0 0 Q_ 

10 0 
2 0 0 

0 0- - - Q 

0 3 0 

0 4 0 

O 0 QL 

0 0 5 

0 0 6 



PHI 

20 

29 30 

31 32 33 

PS1 • 




-A 9- 

13 14 15 16 11 18 19 20 21 



POPULATION 4 



LAMBDA 

0 


0 


0 




I 


0 


0 




2 


0 


0 




0 


0. . 


0 




0 


3 


6 




0 


4 


0 




0 


0 


Ci . _ 




0 


0 


5 




0 


0 


6 




PHI 








34 








35 


36 






37 


36 


39 




PS1 








13 


14 


15 i6 17 18 19 20 21 






3 ' 



-AID- 



INITIAL SOLUTION 
POPULATION 1 



L AMBDA 





1 

2 




1 

0 . 798 
0 . 395 


2 

0.0 
0 . 0 


3 

0.0 

0*0 _ 






- 

3 




0 . 465 


0 . 0 


0 . 0 






4 




0 . 0 


0.796 


0.0 






5 




0 . 0 


0*8 49 


0*0 _ 






6 0* 0 0*75 0 0*0 

7 0*0 0*0 0*597 

R 0* Q 0* 0 0.477 



9 0*0 0.0 0* 538 



PHI 



1 


1 

1.448 


2 


3 












? 

3 


0. 471 
0. 677 


1.127 
0. 249 


1.1 17 












P 51 


1 


2 


3 


4 


5 


6 


7 


fl 9 


1 


0.440 


1.116 


0.872 


0.565 


0.527 


0 . 664 


0.959 


C. 860 0.395 


POPULATION 2 



















LAMBDA 



1 2 3 

1 0. 790 0.0 0.0 

2 0. 395 0. 0 0.0 

3 0* 465 0,0 0*0 

A 0.0 0* 796 0.0 

_5 0.0 O. 849 0.0 

t 0*0 0.750 0.0 

7 0*0 0.0 0*597 

fl 0.0 O* 0 0.477 

o 0* 0 0. 0 0.538 

_ PHI 



1 2 3 

- 1 : 

7 0* 493 I* 03 7 

3 0. 106 Or 208 0.7 17 



PST 

1 2 3 4 _ 5 6 7 6 ? 

‘ 1 0.440 1. 116 0.872 0.565 0.527 0.664 0.959 0. 660 0. 895 

POPUL ATION 3 — 



L AMBDA 



1 2 3 

1 0.798 0.0 0.0 

2 0* 395 0. 0 0* 0 



O 

ERLC 



38 



O' c 



-Ali- 



3 0.465 0,0 0,0 

4 0.0 0.796 0.0 

£ 0.0 0.649 0.0 

6 0.0 O. 75 0 0.0 

7 0.0 0.0 0. 597 

e 0.0 0. O 0.477 

9 0.0 0.0 0.538 

PHI 



12 3 

1 0.956 

2 0.582 0. 933 

3 0. 560 0.485 0.845 



PSI 

1 2 3 4 5 6 7 8 9 

1 0.440 1.116 0.872 0*565 0.527 0.664 0.959 0.860 0.895 

>UL ATIQN 4 „ 



LAMBDA 



1 2 3 

1 0. 798 0. C 0.0 

? 0. 395 0. 0 0^0 

3 0.465 0.0 0.0 

4 0.0 0.796 0.0 

5 0. 0 0. S49 O.Q 

6 0.0 0.750 0.0 

7 0.0 0. 0 0.597 

JJ 0. 0 0. 0 0.477 

9 0.0 0. 0 0.538 

PHI 



12 3 

X 0. 878 

2 0. 483 0.880 

3 0.999 0.363 1.330 



PSI 

1 J> 3 4 5 6 7 8 ? 

1 0.440 1.116 0.872 0.565 . 0.527 ' 0.664 6.959 0.860 0.895 



39 




-A12- 



BEHAVIGR UNDER STEEPEST DESCENT ITERATIONS 



ITER TRY 



ABSCISSA 



SLOPE 



FUNCTION 



0 0.0 

1 O.IOOCOOOOD 00 

2 0, 2466 752 BD 00 



-0.23365917D 03 
-0.13893578D 03 
— O. 493682 80 D 02 



0.IA5A9722P 03 
0. 12720243D 03 
0*1138 B437D 03 



0.39349P76D 00 



O. 116046440 02 



0. 1 11352380 



2 



3 



A 



0 Q. Q ~0t 8 5863 19 AD 02 , 0. 111352380 03 

1 0.393*98760 00 0*21*666360 03 0»12774018D 03 

2 0.135115690 00 0*134910580 01 0*104407170 03 



0 0.0 -0. 555000310 02 0.104407170 03 

1 0. 15511 5690 00 C.38443661D 02 0.103709960 03 

2 fU 816778010-01 -O, 431493360^01 _ 102243 7 SQ__Q^ 



0 

_L 



0-5 n 

0 . 81 6-7.7-ao.lPcQl. 



-0.438528850 02 
-Q* ] 



0. 102243750 03 




40 



-A13- 



BEHAVIOR UNDER FtEPOW ITERATIONS 


ITER 


TRY 


ABSCISSA 


SLOPE 


FUNCTION 




1 


0 


0. 0 


—0. 27' 14 56 08 D 02 


Co 10031534D 


03 




1 


0. lOOCOOOOD 00 


-0*243450030 02 


0.97741812D 


02 




2 


0. 1306951 7D 01 


0. 2 1946 2 6 5D 02 


0.908821990 


02 




3 


0. 892 8471 2D 00 


-0.10321629D 01 


0. 8713331 9D 


02 


7 


0 


t 

o 

• 

° 


-0.44337646D Oi 


0. 8713331 9 D 


02 




1 


0. 892 8471 2D 00 


0. 2 3624045D 02 


1578710D 


02 




2 


0* 3 51 46 57 CD 00 


0*680901460 00 


0. 86359070D 


02 




3 


0. 31668022 D 00 


-0.18325592D-01 


0. 86347 700 D 


02 


3 


0 


0. 0 


— 0. 3785076 ID 01 


0. 86 3477 00D 


02 




1 


0* 3166 8C22D 00 


0.87552759D 01 


0.871724610 


02 




2 


0. 918617420-01 


0. 444271 33D-0 1 


0. 861 782 50 D 


02 


4 


0 


0 

1 

o 


— 0. 1 7366576D 01 


0. 86 1782 50 D 


02 




1 


0* 91861 742 D — 01 


0.15386002D 00 


0. 86 1 0 5439 D 


02 


5 


0 


o 


— 0. 51 584948D 00 


0. 861 05439D 


02 __ 




1 


0. 91861742D-01 


0. 120097280 01 


0. 8613555CD 


02 




2 


0* 286191 890-01 


0. 26 3741 31 D— 03 


0. S6098023D 


02 


6 


0 


0* 0 


-0.37455577D 00 


G, 860900230 


~oT 




1 


0. 286191 B9D-C1 


0.39157427D 00 


0. 860982 93D 


02 




2 


0. 1 393 97 7 00— 01 


0. 1 369632 ID-05 


P * 8609 541 6D 


02 


7 


0 


0. 0 


-0. 23740794D 00 


0. 860954 16 D 


02 






0. 1393977CD-01 


-0. 12891634D— 02 


0. 86093752 D 


02 


8 


0 


0. 0 


-0. 15 794796D 00 


0.8 609 37 52D 


02 




l 


0. 1 393^7700-01 


-0. 38015644D-01 


0. 8609238 7D 


02 




i 


0. 18 37566 ID— 01 


0.38591 097D-05 


0. 86092302D 


02 


9 


3 


0. 0 


-0.21499105D 00 


0. 860923020 


02 




1 


6. 1 6375661D-01 


0. 24735026D-01 


0. 8609 0 55 2 D 


02 




2 


0. 1648578 ID— 01 


— 0.6 977 86 2 2D— 06 


0. 86090 52 8D 


02 


10 


o 


0. 0 


-0.49997638D 00 


0.860905280 


02 




1 


0. 1 648578 ID— 01 


-0.2063L646D 00 


0. 860 84703D 


02 




2 __ 


0- 2 797 93 6 3D — 01 


0.2 7 43 1202D— 04 


0. 86083516P 


OZ_ 


l 1 


0 


0. 0 


-0.726083180 00 


0. 86083516D 


02 




1 


0. 279793630-01 


- 0. 759500480-01 


0, 860 72267 D 


02 




2 


0. 312132 8 7D— G1 


0.1 16 28 02 0D-04 


0. 86072 1440 


02 


12 


0 


0* 0 


-0. 362 40696D 00 


0. 86072 144D 


G2 




i 


0. 3121328 70 — C'i 


0.48471324D 00 


0. 86074044D 


02 




2 


0.1 336797 00 '01 


0.506 31 263D-05 


0. 8606972 1 D 


02 


13 


0 


0. 0 


~~ ~0. 144536 1 4D 00 


0. 8606 97 2 1 D 


02 




i 


0. 133679700-01 


0.2133161 ?D 00 


0. 8607 01 79 D 


02 




2 


0. 54061 C48D-02 


0.4 76 53 9 09 D— 06 


0. 8606 9330 D 


_0_ 2 


14 


0 


0 . 0 


-0. 43866881D-01 


0.860693300 


02 






0. 54061C48D— 02 


0. 1199 40 2 0 D 00 


0. 0 606 9536 D 


02 




2 ~ 


0. 14501964D— 02 


0. 843 5884 8D— 07 


0.860692980 


02 


..15 


0 ._ 


0 . 0 


-0. 15625268D-01 


0. 8606929 8 D 


02 



-Aik- 





1 


Or, 145C1S64D-02 


0. 184174960-01 


0# 86 06 9300 D 


02 




2 


0. 6 65 302 9 40— 03 


0.152483430-09 


0# 860692930 


02 


16 


0 


0. 0 


-0. 132036270-01 


0* 860692930 


02 




1 


0. 665302940-03 


0. 149639240-02 


0. 860692090 


02 




2 


Oe 59757824D— 03 


0. 187024630-07 


0.860692890 


02 


17 


0 


0. 0 


-0. 336 932 71 D— 0 1 


0. 860692890 


02 




1 


0. 597 57 8240—03 


-0. 189336270-01 


0. 860692740 


02 




2 


0. 13641498D-02 


0. 185598170-05 


0. 860692660 


02 


IS 


0 


o 

« 

o 


-0. 299321660-01 


0.860692660 


02 




1 


0. 136414980-02 


-0.66 0297220-03 


0. 860692450 


02 


1 9 


0 


o 

« 

o 


-0. 193316030-01 


0. 860692450 


02 




l 


0. 1364149 80— 02 


0.117416520 00 


0.860693120 


02 




2 


0. 1 S2 8 7C31D-03 


0. 656447340-08 


0. 860692440 


02 


20 


0 


0. 0 


-0. 197791900-02 


0. 860692440 


02 




1 


0. 19287 C3 10-03 


0. 210834990-02 


0. 860692440 


02 




2 


0. 93357C1 CO — 04 


-0. 16 0 824 9 3D— 07 


0. 860692430 


02 


21 


0 


0. 0 


-0. 29334193 0— 0 2 


0. 860692430 


02 




1 


0.9335701 CD-04 


-0. 167599040-03 


0. 860692430 


02 


22 


c 


0. o 


-0. 685567270-02 


0. 360692430 


02 




1 


0. 9335701 CO- 04 


-0.424336540-0^ 


0. 86069243P 02 . 




2 


0. 245003760-03 


-0.222160350-06 


0. 860692420 


02 


23 


0 


o 

« 

o 


-0. 326320370-02 


0. 860692420 


02 




1 


0. 2 4 5CC3 760-03 


0.773092230-02 


0. 860692430 


02 




2 


0. 7 272 C3 930-04 


0. 478375060-07 


0. 860692420 


02 


TIME- ’ 


57. 


37 









4 ^: 




-A15- 

MA X I MUK L I K E L I HO 00 50 L UT I ON 

INP= 0 



POPULATION 1 

Uambda 



1 2 3 





1 

2 

3 


0. 79 8 
C. A? 3 
0. 568 


0. 0 
0. 0 
0. 0 


0.0 

0.0 

0.0 






' -■ " 


A 


0. 0 


0. 796 


0.0 








5 


0. 0 


0. 8A7 


0.0 








6 


0. 0 


0. 7A6 


0.0 




. 


— 


7 


0. 0 


0 . 0 


0.597 








8 


0. 0. 


0 . 0 


0. A96 








9 


0 . 0 


0 . 0 


0.578 


. - 


— 



PHI 



— 




1 


2 


3 






1 

2 




0. 83 8 
0. A 7 5 


0. 911 








~~ 3 

P S I 




0. 908 


0. 348 


1.238 






I 




1 

0. 693 


2 

0. 90A 


3 

0.858 


A 

0 .604 . 


5 6 7 

0.530 0.665 0.812 



POPULATION 2 



lambda 


1 


2 


3 










1 


0. 798 


0. 0 


0.0 










2 


0 . 47 3 


0. 0 


0.0 










3 


0 . 568 


0. 0 


0. 0 ... 










A 


0. 0 


0 . 796 


0.0 










5 


0, 0 


0 . 8 A 7 


0.0 










6 


0. 0 


0. 746 


0.0 







.. 




7 


0. 0 


6. 6 


0.597 










8 


0. 0 


0. 0 


0.496 










9 


0. 0 


0. 0 


0.578 




. 




— 



8 

0. 875 



OH I 



1 '2 3 

1 0.731 

2 0 . 662 C* 92 8 

3 0.522 0.4B6 0.557 



P.SJL 



1 2 3 4*5 6 7 

0. 693 0. 9 04 0, 35 8 O.frOA, Q»53_0 Oj» 66.5 0.JELL2 

POPULATION 3 



a 

0, 875 



LAMBDA 



1 2 3 




4 u 



9 

0.825 



9 

0.825 



^16- 



1 


0.798 


0. 0 


0.0 








2 


0.473 


0. 0 


0. 0 










O- 56 8 


0. 0 


0.0 








- 2 . 

4 


o.o 


0. 796 


0. 0 








5 


c.o 


0. 847 


0. 0 










0* 0 


0. 746 


0.0 






. ... 




7 


V- 

0. 0 


0. 0 


0. 597 








8 


0. 0 


0. 0 


0.496 








9 


0. 0 


0. 0 


0. 578 




— - 





phi 





1 


2 


3 












1 


0. 592 
















2 


0. 469 


1. C69 










„ 


- - - -- 


3 


0.074 


0. 199 


0. 853 












p SJ 




















1 


2 


3 


4 


5 


6 


7 


8 


1 


0. 693 


0. 904 


0. 858 


0.604 


0.530 


0.665 


0.812 


0*8 7 5 


ON 4 


















LAMBDA 




















1 


2 


3 












1 


0.79 8 


0. 0 


0. 0 












2 


0.473 


0. 0 


0.0 












3 


0, 568 


0. 0 


0. 0 










... 


4 


0. 0 


0. 7S6 


0.0 












5 


0. 0 


0. 847 


0. 0 












6 


o. o 


0. 746 


0.0 










... . 


7 " 


0.0 


>0* 0 


0. 5 97 












8 


0. 0 


0. 0 


0.496 












9 


0. 0 _ . 


0 . 0 


O. 578 








— 


. 



PHI 



1 

2 


l 

1. 12 8 
0. 386 


2 

1* 127 


3 












3 

ps r 


0. 6 27 


0. 265 


1.191 














1 

0. 693 


2 

0. 904 


3 

0. 858 


4 

0.604 


5 

U. 530 


6 

0«» 665 


7 

0.812 


8 

0.875 



9 

0.825 



4/1 




O' O; 



■AIT 



TEST HF GOODNESS OF FIT 

CHI SQUARE Wl TH 141 DEGREES OF FREEDOM fs 172- 1385^ 

PROBABILITY LEVEL IS C.Q38 




4o 



-4 C7 4 ui t-P* uj r\> V- 



POPULATION 1 



S.I. GMA = LAMBDA *_P HI » LAMBDA *+PSI ** 2 



1 2 3 



0. 3 16 " 1. 004 

0. 379 0. 225 1.006 

0. 302 0 , 179 0.215 

0.321 0. 1 91 0. 2 29 

0. 283 0. 168 0.201 

0.^33 0. 257 0. 30 8 

0. 359 0. 213 0.256 

0.419 0. 2 49 0.298 



RESIDUALS = SIGMA-S 



4 



5 



0.942 

0.614 0.934 

0.541 0.575 

-0 * 16 6 0.17 6 

0.137 0.146 

0.160 0.171 



6 



7 



0.949 

0_-155 .1.10' 

0.129 0.36; 

0.150 0*428 



8 



1.070 

0.355 



9 



1 - 



— 


— 


1 


. _ _2 3 


4 


5 


A 






1 

2 

3 


-6. 053 
0. 031 
0. 007 


0. 2 59 

0.JL6 0 -0.120 






O 


_ _7- 




4 

5 

6 


-0. 021 
0. 11 3 
-0. CO 9 


-C.030 -0.117 

0. 064 —0. 054 

- _ -0. 013 -0. 155 


-0.076 
0.040 
—0. 06 1 


0.090 

— ri n it 








7 

8 
_9_ 


- 0. 056 
0. 143 
-0. 044 


-0. 007 -6. 027 

0* 204 0. 089 

-- -0. 009 -0. 155 


-0.014 
0.074 
0.04 2 


VJ • U 1 1 

-0.005 

-0.036 

0.102 


-0.1 48 
-0. 145 
-0.035 


0. 128 
-0.005 














— 0 • 0 51 


0.059 


































































































































— ■ 


- . 


... 













8 



9 



- 0.022 

-0.136 -0.- 




46 



-A1Q- 

POPULATTOM 2 





SIGMA = LAMBDA^PHI *1 AMBDA 1 -i-PSI **2 
















1 


2 


3 


4 


5 


6 


7 


8 


9 


1 


0, <345 


















2 


0* 2 76 


0*981 
















3 


0. 331 


0- 196 


0. 972 














4 


0. 357 


0-212 


0-2 54 


0-953 












5 ‘ 


0. 380 


0- 226 


0-270 


0-626 


0-946 










6 


0. 334 


0- 199 


0.238 


0-551 


0.586 


0- 958 








7 


0. 249 


0- 148 


0. 1 77 


0-2 3.1 


0.246 


0-2_17 


0.853 






p 


0.206 


0.122 


0. 147 


0- 192 


0.204 


0* 180 


0. 165 


0.903 




Q 


0. 241 


0-143 


0- 1 71 


0-224 


0-238 


0*210 


0. L92 


0. 160 


0. 867 




HP SlriuAL5~=“SIGMA— S 




















1 


2 


3 


4 


5 


6 


7 


8 


9 


1 


0. 049 


















2 


-0.C57 


-0- 093 
















3 


- Q.J026 


-0. 004 


0. 126„ . 














4 


0. 001 


-0* 117 


-0* 029 


-c70 30 












5 


0. 017 


0- C6 7 


0- 059 


•0-030 


0.026 










6 


- 0. 044 


0. 070 


— 0_. 06 9 


0-032 


-0.007 


0. 052 








7 


-0. 021 


-0. 080 


-0. 048 


0-094 


0 • 090 


-0.023 


0.197 






8 


C. C80 


0-114 


0* 078 


-0 -064 


-0.035 


-0.089 


-0.045 


0. 211 




9 


r P. 0*^3 _ 


-0. 105 


__ Q.098 „ 


-0.075 


0*009 


-0*018 


-0 081 


-0. 120 


0. 083 




47 



to O' ui ro »- ; o oc ,*4 o j$* v» r\> * 



lawb oa^ phi ^lambda^psx ** g 



SI GMA ,._= 



I 


? 


3 


4 


5 


6 


7 


0. 856 














0. 224 


0. 949 












0. 268 


0. 159 


0. 927 










C„ 298 


0. 1 77 


0. 212 


1.036 








0. 317 


0. 1 88 


0. 225 


0.714 


1.041 






C. 279 


0. 166 


0. 198 


0.629 


0.669 


1 c 0 31 




0. 035 


.... .P.021 


0„ 025 


0.095 


0. 101 


0,089 


0, 964 


0. 029 


0. 017 


0.021 


0.079 


0.084 


0.074 


0. 253 


0. 0 34 


0. 020 


0.024 


0.092 


0.097 


0.066 


0. 295 


OUALS = 5 1 GMA— S 












1. 


2 


3 


4 


5 


6 


7 


-0. 066 














0. 025 


0.204 












__ -D*0A9 


— 0 . 029 


- 0.^0 5 4 __ 










6. ooi 


0. 1 35 


-0. 003 


o7io¥~ 








-0. 040 


-0. 023 


0. 036 


0.018 


-0.086 






~Pjl£70_ 


0. 07ft 


Qj^0_8Q .. . . 


0.054 


„ .-0, 038. 


0,013 




-0. 179 


(J. 012 


0. 095 


0.007 


-0.018 


—0. 034 


-0.063 


0. 051 


0. 027 


0. 164 


0.025 


-0.011 


0. D 40 


0.039 


-0,_052_.._ 


0. 1 41 _ 


0 ..084 


-0.063 


0.076 


— 0. 035 


0. 143 




— 





— 


— — — 








8 



0.975 

0.245 



e 



-0. 254 
-0.077 



9 



0. 966 



9 



-0. 032 
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POPULATION 4 



-A21- 





SIGMA - LAMBDA*PHI*LAMBDA ■ 


*-PSI**2 














1 


2 


3 


4 


5 


6 


7 


8 . 


1 


1*197 
















2 


0. A 2 6 


1.070 














3 


0. 51 1 


0.303 


1 . too 












4 


0. 245 


0. 146 


0.175 


1.079 










5 


0. 26 1 


0. 155 


0.186 


0.760 


1.089 








6 


0. 230 


0. 136 


0. 163 


0.669 


0.712 


1.069 






7 


0. 299 


0. 177 


0.2 13 


0.126 


0.134 


0..1 18_ 


1.085 




8 


0. 248 


0. 147 


0.176 


0.105 


0. Ill 


0.098 


0.353 


1. 058 


9 


0. 289 


0. 172 


0.206 


0.122 


0. 130 


0.1 14 


0.411 


0.341 




RESIDUALS = SIGMA-S 




■ 


— ” - • - 


. 










1 


2 


3 


4 


5 


6 


7 


8 


1 


0. 072 
















2 


0. 016 


-0. 391 














3 


-0.012 


-0. 106 


0. 047 












* 


-0. 062 


0. 1 33 


0. 1 1 1 


oToo 7 










5 


-0* 032 


0. 14? 


0. 175 


-0.064 


-0.038 








6 


-0. 190 


-0. 175 


0. 062 


0.054 


_0_. 050 


0.089 






7 


-0. 224 


-0. 277 


-0. 052 


-0.056 


0.047 


-0*300 


-0.292 




01 


0. 120 


0. C67 


0. 145 


0. 183 


0.047 


-0.091 


0.012 


. 1’54 


9 


0.034 


0.222 


0. 195 


_Q.jJL7.6_..- 


0.019 


-0* 134 


0. 178 


.05 7 




- — - 


— 





• - 







- 
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-A22 



SCALED SOLUTION! 
PO PUL AT I ON 1 

L PDA 





1 


2 


3 




1 


0. 721 


0. 0 


0.0 




2 


0.428 


0. 0 


0. 0 




3 


0. 513 


0. 0 


0.0 




/• 


0. 0 


0.797 


0.0 






D. 0 


0. 848 


o.c 




6 


D.O 


0. 74 7 


0.0 




7 


D. 0 


0. 0 


0. 583 




6 


0. 0 


0. 0 


0. 484 




9 


0. D 


0. 0 


0.5 65 





PM [ 



1 ? 3 

.i _..i*.a24 

2 0. 52 0. 90P 

3 l.C 0.356 1.298 



P SI 



1 



i 

O. 6*33 



0. 9 04 



3 

0. 858 



4 

0.6D4 



5 

0.5 30 



(Sl . 

0.665 



7 

0 © 8 1 2 



8 

0.875 



9 

O. 82 8 



POPULATION . 2 



LAMP DA 







1 


7 


3 




1 




0. 7 ? 1 


0. 0 


D. 0 




2 




0. 4?8 


o._o 


0.0 




3 




0. 513 


0. 0 


0. 0 




4 




0. 0 


0. 7C7 


D. 0 




5 




0*0 


Or. 84 8 


D.D 




6 




o. 6 


0. 74 7 


D.O 




7 




0. 0 


0. 0 


0. 583 




9 




0. 0 


0. 0 


0. 434 




o 




0. 0 


0. 0 


0. 565 





PHI 



1 2 3 

1 ,_.a« d.9fi 

2 0. 621 0. 975 

3 C. 591 0.497 0.584 



PS I 



1 

POPUL Al l ON 3 



1 2 3 4 5 6 7 

0.693 0.904 " 0.858 0.604 0.530 0.665 0.612 



8 

0. 875 



9 

0. 825 



L A V 8 Q A 

1 

1 0.721 

2 0*4 2 H 



2 ' 3 

0.0 O.D 

_0*_Q 0.0 



O 

ERJC 



-A 23 - 



3 

4 

_5 

6 

7 

it 

<5 



0.513 0.0 0.0 

0.0 0. 7 97 0.0 

o. 0 Q»D. 

0.0 0.747 0.0 

C. 0 0 . 0 0. 583 

0. 0 0.0 0.464 

0 . 0 0 . 0 0 6 5 



PH I 



2 

3 



12 ’ 

il..724__ „ 

0.518 1.056 

0. C 84 0.203 0.895 



PSI 



\ 2 3 4 5 6 7 8 

1 0.693 0. 904 0.858 0.604 0.530 0.665 0.512 0.875 

L ATIQM 4 





L &M 6 DA 












1 


2 


3 


1 




0. 721 


0. C 


0.0 


2 




0. 4 25 


0. 0 


0. 0 


3 




0. 513 


0. c 


0.0 


4 




0. 0 


0. 797 


0. 0 


5 




0. 0 


0. 84_9_ 


0. 0 


6 




0. 0 


0. 7 3 7 


G.O** 


7 




0. 0 


0. 0 


0. 583 


p 




c. 0 


_0j, 0 


_0*_4_84 


9 


PHI 


c. 0 


0. c 


0.565 






1 


2 


3 


. ,v_ 




1. 380 






2 




0. 426 


1.124 




3 




0. 7 1 0 


0. 271 


1.249 



DSI 



1., 2 3 4 _ 5 6 7 8 

1 0.^93 0.9 04* 0 • 8 5 8 *0.604 0.530 0.665 0.812 0.875 



bi 




9 

0. 825 



o> 6 



- A 24 - 



STANDARD ERRORS 



P O PUL AT ION 1 



LAMBDA 



12 3 

1 0 . 0 0.0 0.0 

2 0. C83 0. 0 0.0 

3 0.087 0.0 0.0 

A- 0.0 O.C 0.0 

5 0. 0 O o 06 2 0. 0 

6 0*0 0.059 0.0 

7 0.0 0.0 O.C 

B, 0_#,Q OjO O. 091 

9 0.0 O.C 0.098 

PHI . 



1 2 3 

. 1 0 „ 2 A ? 

2 0. 1 5? 0. 195 

3 0.232 0.177 0.396 



PS I 

1 .? 3 4 5 6 7 

1 ~ 0.055 '.041 0.042 " 0*037 " 0.042 0.036 0.047 

P OPUL, AT I ON 2 

lambda 





1 


p 


3 


l 


0. 0 


0. 0 


0. 0 


p 


__ _0a,C : 3 ... 


a, o 


0.0 


3 


0. C 8 7 


0. 0 


0.0 


4 


0. 0 


0. 0 


0.0 


5 


0. 0 _ . 


_0. 062 


0. Q 


6 


0.0 


0. 059 


0-, o 


7 


0. o 


0. 0 


0. 0 


R _ 


0. 0 


0. 0 


0. 091 


9 


0. 0 


o 

! • 

IO 


0. 098 



PH I 



1 2 3 

. 1 . - 0*^.2 0 _ 

2 0.151 0.196 

3 0.163 0.154 0.237 



s i 

_1_ 2 _J 4 _ 5 6 

1 ~ 0.055 0.041 0.042 ~ 0.037 6.042 0.036 

POPULATION .^ 



LAMBDA 



12 3 

1 0.0 0.0 0.0 

2 _0.C8 3 0._0 3^0 



7 - _ 

0.047 




52 



8 

O. 043 



8 

0. 043 



O j ) o & 



-A25- 



3 0. 037 0* C 0*0 

A > 0. O 0*0 0,0 

0*0 0*062 0*0 

6 0*0 0 * 059 0.0 

7 0.0 0*0 0*0 

jj Q.Q _Q* -0 0* 09 .1 

g 0.0 ^T5 0.098 



PHI. 



_ JL - 

? 






1 2 3 

J^ lZDQ - ~ 

O ^ H '? 0.226 

C 5 1 5 A 0. 167 0.311 



psT 



l 



1 2 3 A 5 

0*055 0.041 0.042 0*037 0*042 



6 7 

0.0~ b 0.047 



UL A T IQ N A 



L AMBDA 






1 


2 


3 


1 


0. 0 


0. 0 


n. o 


2 


0* CP 3. 


0. 0 


.0.0 - 




0. C 8 * 


0.0 


o • 6 


4 


0* 0 


0. 0 


0.0 


Cj 


_._0 % o 


_Q, 062 


0.0. .. _ _ _ ..... . . _ » — — -■ - • 


6 


0* 0 


0.05<= 


o.o 


7 


0. 0 


0. 0 


0.0 






r,. p 


0 * 09 X - - - • 


9 




0.0 


o.o 


0* 09R 


P MI 






_ - . 




1 


2 


3 




0. 







" 2 
3 


0. 182 
0. 233 


6.2 43 
0. 194 


0.395 



P ST 



1 Z 3 A 5 6 7. 

67 0$5 0. 041 0 * 04 2 0*0 37 0*042 0*0 36 0.04 7 



8 

0.043 



0. 043 





9 

0. 046 



O J> 



■A26 



S I MULTAN 10 US FACTOR ANALYSIS IN SEVERAL POPULATIONS 

ART I FI CtAL DATA FOR ILLUSTRATIVE PURPOSES 

NIP ( 1 >~ 61 

~Nr>( 2T^ isa" 

2 

1 0 
3 

LOGI CAL IIVCICATOaS t OOLUMNS 5 1-57) sTTFTF FF 

OUTPUT I NO IC A '“'1^5— l l 

ESTIMATED TIME IN SECONDS^ 220* 



5 ^ 




-A2.J- 



UL AT \ OH I 



1 

z 


1 

i. ce 8 

0. 254 


2 

1.124 


3 


4 




6 


7 


"a ' ' 9 


3 

4 

c t 


-0. 050 

0 . oil 

^C. 0 54 


0. 408 
0. 522 
0.307 


1.448 
0.4 07 
0. 012 


1 n 05 2 
0.064 


1.071 






1.3 58 

0.337 0.997 


6 

T 

_ 3 _ 


0. 1 1 0 
0.248 
0.231 


0. 291 
0.420 
0. 519 


0. 013 
0.310 
0. 449 


0.011 

0.102 

0.263 


0.820 

0.615 

0.181 


1. 116 
0.659 

0_. 0 86 

n.. 063 


0.981 

0-415 

0.188 


9 

UL AT I ON 2 


0. 396 


0. 127 


Oi 060 


0.031 


-0.083' 


5 
















1 


2 


3 


4 


5 


6 


7 


8 9 


1 


0. 996 
















z 


.0_._P_86^ 


0. 921 














3 


-0. 1 20 


0. 198 


0. 739 






■ ■ 


'■ 




4 


-0. 059 


0. 219 


0.187 


0.981 










5. 


_ _ 0 . 1 54 


0. 2 96 


Oj 041 


0 _ P Z 23 ^ 


0.931 








6 


0. 021 


0.355 


0. 2 09 


0.188 


0.693 


i \ : i6 






7 


0.121 


0. 34 9 


0. 087 


0.110 


0.575 


0. 704 


1. 019 




8 


3* 1 51 


0.212 „ 


0 . 009 


- 0.070 


0.087 


0 . H 7 


0.122 


1 . 01 3 


9 


G. 320 


- 0. 021 


- 0 . 009 


- 0.142 


0.053 


0 . C 93 


0.033 


0.211 1.221 








______ . 
















bo 



o 

ERIC 



-A 28 - 



PARAMETER 


SPECIFICAT IONS 




POPULATION 1 


L AMBDA 




0 


0 


1 


0 


0 


0 


2 


0 


0 


3 


0 


0 


0 


0 


0 


0 


4 


0 


0 


5 


0 


0 


0 


0 


0 


0 


6 


PHI 






7 






8 


Q 




10 


11 


12 


PS I 






1 3 


14 


15 16 17 18 19 20 21 




POPULATION 2 


L AM 8DA 




0 


0 


22 


0 


0 


0 


2 


0 

p__ 


0 

0 


0 


c 


0 


0 


4 


0 


0 


5 


0 


0 


0 


0 


0 


0 


23 


P H I 






24 






2 5 


2*5 





27 28 2<? 



PSl 

30 31 32 33 34 35 36 37 38 

T I ME— 58, 91 



* 







5b 



-A 2Q- 



MAXIMUM L I KB L 1 HO CD SOLUTION 

INP - Q 

’UI.&TiON 1 

I A"H80A 



1 i .... C 

2 < ’OC 

3 ■ 

4 O* 504 

5 0- C- 

0 . 0 

7 0 . 0 

ft 0-0 

9 Q *_0 



2 3 



0.0 


0.485 


0. c 


0.0 


0* 0 


0© 0 


0. ' 


0 © 0 


1 • 000 


0 . 0 


lj 174 


0 . 0 


o . 9 1 : 


0-0 


0- 0 


1 - OOt . 


0*0, 


. 0L*_5 10_ 



PHI 



1 2 3 

1 0-902 

.2 O- 276 0- 64 7 

3 0-473 0-159 0-669 



PS I 



12345678 
1 0, 965 0# 466 1-110 0-870 j0-_5 8 4 0.49 6 0.670 0-830 



?UL ATI ON ? 



LAMBDA 



1 2 3 

1 C-0 0-0 1-527 

2 1 - 000 0-0 0-0 

.3 0- 49 5 0-0 GUO 

4 0- 584 0-0 Or 0 

5 0. 0 1 - 000 0.0 

6 0. 0 1-174 0«_£ 

7 0- 0 0.962 0-0 

3 C-0 0-0 1-000 

9 _0. 0. Q- 0 1 . 984 



P H I 



12 3 

3 0- 408 

2 0.* 30 0_ 0*_6 04 

3 0-004 0-047 0-105 



P SI 



1 2 3 4 5 6 7 8 

1 _ _ 0 9.6-7— 0- 7JL9 0, 7 99 0-91 5 0.J596 0-526 0.-6 62 0.953 



57 




o & j 1 o 



-AT- 
TEST OF GOODNESS OF F I T ~ “ ~~ ~ 

CHI SQUARE WITH 52 DEGREeS OF FREEDOM* IS 88.7 335 

PROBABILITY LEVEL IS O.OOl 1 ... 



A 



f 



o 

ERIC 



5 6 



